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We prove a Carleman estimate for elliptic second order partial differential 
expressions with Lipschitz continuous coefficients. The Carleman estimate is 
valid for any complex-valued function u E W 2,2 with support in a punctured ball 
of arbitrary radius. The novelty of this Carleman estimate is that we establish 
an explicit dependence to the Lipschitz and ellipticity constants, the dimension 
of the space and the radius of the ball. In particular we provide a uniform and 
quantitative bound on the weight function for a class of elliptic operators given 
explicitly in terms of ellipticity and Lipschitz constant. 


1. Introduction 


Carleman estimates were first introduced by Carleman [Car39] in 1939 in order to establish a 
unique continuation property for elliptic operators L with non-analytic coefficients, i.e. if the 
solution u of Lu = 0 in C vanishes in a non-empty open set, then u is identically zero. 
While Carleman’s original result applies to the case d = 2 and L = —A + V with V E L[^ C (M 2 ), 
by now there is a wealth of results pertaining to Carleman estimates and its application 
to unique continuation, see, e.g., [Hor89]. In particular, notable attention has been paid 
to the case V E L[^ c (M rf ), since this can be used to show that the Schrodinger operator 
H = —A + V with potential V E L^ oc (R d ) has no positive eigenvalues, see [JK85] and the 
references therein. Further applications of Carleman estimates are, for example, uniqueness 
properties of solutions of Schrodinger equations [KPV03, IK06, EKPV12], uniqueness and 
stability of inverse problems [Klil3, Klil5], or control theory of partial differential equations 
[FI96, CZ01, Roul2, RL12], 

So far, all the mentioned applications of Carleman estimates concern qualitative statements. 
For its proofs, the particular dependence on the parameters is not essential. Certainly, there 
are plenty of applications which require some knowledge on the dependence on the parameters 
entering in the Carleman estimate. For this purpose, we cite the following Carleman estimate 
by Escauriaza and Vessella [EV03], see also [KSU11]. Let L = — div(MV) be elliptic with 
Lipschitz continuous coefficients and denote the corresponding constants by f}\ and i?2j 
respectively. Then there exist k E (0,1) and C, H, ao > 1, depending on di and i? 2 , and 
a function w : — > R_|_ satisfying |x|/S < w(x) < £|a;|, such that for all real-valued 

u E C£°(B k \ {0}) and all a > «o we have 


a 


w 1 - 2 a |Vii | 2 + a 3 


w - l ~ 2a u 2 <C J w 2 ~ 2a (Lu) 2 . 


(1) 
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Here, B K denotes the open ball with radius k and center zero. 

The exceptional feature of this Carleman estimate is that the second term on the left hand 
side goes with a 3 and the particular powers of the weight function. Thanks to this, Bourgain 
and Kenig [BK05] were able to apply (a refined version of) this Carleman estimate to prove 
that if A u = Vu in M rf , it(0) = 1, |tt| < C and \V\ < C , then for all x £ with \x\ > 1 


max \u(y)\ > c • exp (— c / (log|x|)|x| 4//3 ) . (2) 

Let us note that the quantitative unique continuation principle in Ineq. (2) was crucial 
for the answer to a long-standing problem in the theory of random Schrodinger operators, 
namely Anderson localization for the continuum Anderson model with Bernoulli-distributed 
coupling constants. Moreover, if the precise decay rate 4/3 in Ineq. (2), which results from the 
particular form of Ineq. (1), would be replaced by 1.35, one could not conclude localization 
for the continuum Anderson-Bernoulli model using the same techniques, cf. [BK05, p. 412], 
Ineq. (2) was also used to prove bounds on the density of states measure for Schrodinger 
operators in dimension d £ {1,2,3}, a restriction which stems from the specific parameters 
in Ineq. (2). 

Another example where the Carleman estimate (1) proved to be useful is so-called scale-free 
and quantitative unique continuation. In [RMV13], the authors prove that if |Au| < \ Vu\ in 
Al = (— L/2, L/2) d , then for all 6 £ (0,1/2) one has 

U Wl 2 (A l ) - IMI L 2 (W S (L))’ ( 3 ) 

where C = C(d) and Ws(L ) is some union of equidistributed d-balls in A^. The notion 
scale-free and quantitative refers to the fact that C is independent of L and the dependence 
on 6 is known to be polynomial. The latter fact holds due to the bounds on the weight 
function given above and is essential for their application to random Schrodinger operators. 

The mentioned applications so far concerned Schrodinger operators only. In this case, 
the (abstract) constants k and E in Ineq. (1) depend only on the dimension. Due to the 
method of proof in [BK05, BK13, RMV13], a precise knowledge of k and S would change 
only dimension-dependent constants like c, c' and C. However, if one, for example, pursues to 
prove an analogue of Ineq. (3) for solutions of \Lu\ < \Vu\ with L as above, one encounters 
dependencies of the constants k and S on the ellipticity and Lipschitz parameters of L. If we 
have no quantitative bounds on these dependencies, the use of Carleman estimate (1) might 
not be feasible, cf. Remark 2.4 for details. Such problems have already been discussed in the 
recent proceedings [BTV14] and [BNR + 15] on quantitative scale-free unique continuation 
principles for elliptic operators. This strongly suggests that a refinement of the Carleman 
estimate (1) with explicit bounds on k and E is a worthwhile goal. 

In this note we give a precise refinement of the Carleman estimate (1) in the sense that 

(i) the estimate is valid on the whole punctured unit ball (i.e. k = 1), and 

(ii) all the constants, including the bound E on the weight function, are explicitly calculated 
in terms of ellipticity and Lipschitz constant. 
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Our proof is based on techniques developed in [EV03, BK05]. Compared to the Carleman 
estimates of [EV03] and [BK05], our result holds with k = 1 and additional gradient term in 
the lower bound at the same time. Let us emphasize that, beside the quantitative control 
on all the parameters, even this involves non-trivial modifications of the existing proofs, cf. 
Section 4. 

The usefulness of the established quantitative bounds has been shown by the proof of an 
analogous result to Ineq. (3) for an explicitly given class of elliptic operators, see [BTV15] 
and Remark 2.4. 

The paper is organized as follows. In Section 2 we state the main result of the paper. In 
Section 3 we provide preparatory estimates and identities needed for the proof of the main 
result. In Section 4 we give the proof of the main result, while the proofs of two technical 
lemmata are postponed to the appendix. 


2. Main result 


Let d £ N and L be the second order partial differential expression 

d d 

Lu := — div(AVu) + fe T Vu + cu = — di (a^dju) + bidiU + cu, (4) 

i,j =i i,j =i 

acting on complex-valued functions on M d , where A : —y with A = {oA) d j = \, 

b : JR 6 * — > C d , c : — > C, and di denotes the z-tli weak derivative. Moreover, we denote by 
B p C the open ball in with radius p > 0 and center zero, by \z\ the Euclidean norm 
of z £ C d , and by HMUqo, ||M||i and \\M\\ the row sum, column sum and spectral norm 
of a matrix M £ C dxd . For the coefficient functions A, b and c we introduce the following 
assumption. 

Assumption. Let p > 0, > 1 and d-z > 0. We say that A(p, di,i? 2 ) is satisfied if and 

only if b, c £ L°°(B p ), a lJ = a Jl for all i, j £ {1,..., d} and for almost all x,y £ B p and all 
£ £ R. d we have 

^r 1 |e| 2 <e T ^)e<^i|?| 2 and \\A(x)-A(y)\\ 00 <'& 2 \x-y\. (5) 

By Rademacher’s theorem [Fed96], if A(p,$ i,i? 2 ) is satisfied, then the coefficients a 13 are 
differentiable almost everywhere on B p and the absolute value of the derivative is bounded 
by i? 2 - Operators L for which the lower bound of the first inequality in (5) is satisfied are 
called elliptic. Elliptic means therefore, that for almost every point x £ B p the symmetric 
d x d matrix A is positive definite. A simple example is A = I. In this case the operator L 
coincides with the negative Laplacian. We use the notation Aq = A(0). 

For /i, p > 0 we introduce the function w p . p : —>• [0, oo) by 


w p , p {:x) := (p(a(x/p)), 

where a : W 1 —> [0, oo) and ip : [0, oo) —> [0, oo) are given by 
r(x) := (x t Aq 1 x ') 1//2 and ip(r) := rexp ( — 
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The function ip obeys the upper bound <p(r) < r. For a lower bound we distinguish two cases. 
If \fid\p < 1 we use the lower bound <p(r) > re~^ r . If \/M ip > 1 distinguish two subcases. 
If r < 1/p we use 1 — e^* < pt. and obtain <p(r) > re~^ r > rie^fidfp)^ 1 . If r > l//x we split 
the integral according to the two components [0,1 / p] and [1/p, ?’]. On the first component 
we use 1 — e _/xt < pt, on the second one 1 — e _/it < 1. This way we obtain ip(r) > l/(ep) if 
r > 1///. Putting everything together we obtain for all r G [0, \/Mi] 


<p(r) > —, where p\ = 
Ml 

For the weight function w P)P there follows 

(j[x) 


e'/®!^ if y/W\ p < 1, 
e\fdfp if yfdifj, > 1. 


Mx G B 


p ■ 


PM i 


< ^ < dh < 

PMi P P 


Keep in mind that we will drop the index of the weight function in Section 3 and onwards, 
and write w instead of w PiP . Our main result is the following Carleman estimate. 

Theorem 2.1. Let p> 0, Mi > 1, M 2 > 0, Assumption A(p, Mi, M 2 ) be satisfied and 


p > 33dM^ i// "'$2p. 


Then there exist constants ao = ao(d, p, Mi, M2, m> ||&||oo> ||c||oo); C = C{d, Mi, pi?2, m) > 0, 
suc/i that for all a > ao and all u G VF 2,2 (M rf ) with support in B p \ {0} u>e have 

f (ap 2 w 1 p fj 2a X7u T AX7u + a 3 wf jfi 2 a \u\ 2 ) dx < Cp 4 f iu 2 ^ 2 "|Lu| 2 dx. (6) 

’ ’ JR d 

The constants ao and C are given in Eq. (28). 

Remark 2.2. In the case where b and c are identically zero, the conclusion of Theorem 2.1 
holds with C = C and cto = do where C and do are given in Eq. (26). A straightforward 
calculation shows that they obey the upper bounds 

C < 2d 2 d\e A ^p\ (3/r 2 + (9pM 2 + 3)p + l) Cfi 1 


and 


d 0 < llrf 4 Mf /2 e 6 ^V?(3pM 2 + m + l) 2 (1 + M(M + 

11 /2 

where C p = p — 33' M 2 p. In the general case where 6, c G L°°(B p ) the conclusion of the 
theorem holds with 

(7 = 6(7 and ao = max jdo, C'p 2 |[6||^ 0 M^ 2 , C' 1 / 3 p 4 / 3 ||c||^ 3 i/Mij , 

cf. Eq. (28). We observe that if M 2 0, then both upper bounds of C and ao tend to infinity 
if p tends to 33 dd 1 M 2 P- If $2 = 0 the upper bound of do is uniformly bounded for all 
p > 0. In the case of the pure negative Laplacian, i.e. A = Id and b, c = 0, we have Mi = 1 
and M 2 = 0. In this case we note that our Carleman estimate is valid for arbitrary p > 0. 
For example, in the case p = 1 we infer from the proof of Theorem 2.1 (not using the above 
estimates) that C < 8e 8 M 2 and ao < 18e 12 d 4 . 
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Remark 2.3. Our result gives rise to a uniform Carleman estimate for a certain class of elliptic 
operators given precisely in terms of ellipticity and Lipschitz constants. More precisely, fix 
> 1 and t? 2 , b, c > 0, p > 0 and denote by C the class of elliptic operators of the form (4) 
for which 

(i) the ellipticity and Lipschitz constant on B p are bounded by d\ and $ 2 , and 

(ii) the coefficient functions b and c satisfy || 6 ||oo < b and ||c||oo < c. 

Then, by monotonicity we can apply our result uniformly for all L E £ with, e.g., p = 
34 dd 1 7 P- In particular, there are constants do and C, depending only on d, di, b, c 
and p, and a function w = w P}P : —> [ 0 , oo), such that ( 6 ) holds for all elliptic operators in 
Led 

Remark 2.4. As mentioned in the introduction, the usefulness of our result has already been 
shown in the recent preprint [BTV15]. More precisely, based on the Carleman estimate 
from Theorem 2.1, they show so-called quantitative unique continuation principles, sampling 
theorems as well as equidistribution theorems for solutions of certain variable coefficient 
elliptic partial differential equations or inequalities. In order to stress the usefulness of 
our result, we formulate exemplary a special case of the sampling theorem from [BTV15]. 
Thereafter, we discuss why, e.g., the Carleman estimate from Ineq. (1) is not satisfactory to 
obtain such a result. Beforehand, we introduce some notation. 

Let 5 > 0. We say that a sequence Zj E R rf , j E Z rf , is 5-equidistributed, if 

VjEZ d : B{zj,S) C (-1/2, l/2) d + j. 

Corresponding to a (5-equidistributed sequence Zj E M rf , j E Z rf , we define the set 

5 5 = J B( Zj ,S) cR d . 

j&7, d 


Note that the set S$ depends on the choice of the (5-equidistributed sequence. 

Theorem 2.5 ([BTV15]). Assume 

ei := 1 - 33 d{Vd + 2 )^ 1 / 2 (e ^ /2 + 1 ) 0 2 > 0 . (7) 

Then for all measurable and bounded V : —> M, all ip £ IT 2 , 2 (M d ) and (p E L 2 (M rf ) satisfying 
\Lip\ ^ IV’V’I + ICI almost everywhere on all S E (0,1/2) and all 5-equidistributed sequences 

we have 

U\\s s + IICIIk.^ ^ c sfUc||V ; lliRd) 

where 

f 5 \ ( 1 “ ln£l +ITIlM 3 +G ,2 ||6||^ 0 +||c||^ 3 ) 

c s njC = dil—J (8) 

and d\, and d 3 are constants depending only on di, d 2 , and the dimension d. 
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The constants d±, d 2 and d 3 are given explicitly in [BTV15]. Note that the exponent in 
( 8 ) can be interpreted as an estimate on the vanishing order as studied in [DF 88 , Kuk98]. 
The obtained estimate of vanishing order and its dependence on ellipticity and Lipschitz 
parameters originates from the particular form of our Carleman estimate, including C, ao, 
w, and the cubic behaviour of a. 

Certainly, Theorem 2.5 can be proven by using the Carleman estimate (1) of [EV03]. 
However, then Assumption (7) would involve some function of k and E. Since n and S itself 
are abstract functions of and d 2 , it would be not clear if the corresponding Assumption (7) 
is ever satisfied. Moreover, in [EV03] no monotonicity condition on k and S is given. For 
this reason, it is even not possible to verify the corresponding Assumption (7) by choosing 
$1 and $2 sufficiently small. 

3. Preliminary relations and quantitative estimates 

In this section we provide some preparatory estimates and identities for the proof of Theo¬ 
rem 2.1. In particular, we will prove the quantitative estimates needed for our quantitative 
version of a Carleman estimate. First we introduce some notation. Let Lq be the part of L 
containing only the second order term, i.e. 


L 0 u := — div(AVu). 


For k G {l,...,d} let G M. d be the k-th unit vector. For g G {a, w} and real-valued 

/ G C£°(M d ) we define the functions : R. d —> M, Df : —> M and —>■ given by 



gLpg _ 
Vg^AVg ’ 


D f :=w 


Vf T AVw 

X7w T AX7w 




A 
w > 


and 


,a , = gAVg 

9 ' Vg T AVg' 


Note that F^° refers to the function F^° = — g div(AoV g) / (V g T AqV g) — 1. Moreover, we 
introduce the matrix-valued function —y M. dxd , 

:= Aft A + 2 div(^ oA)-\ (. AD(hf) + [AD(hf)f) . 


Here hf A denotes entrywise multiplication, i.e. hf A := (hf 7 ' 9 )fj =1 , while the divergence 
of a matrix denotes entrywise divergence. The matrix D(h^) is defined as the d x d matrix 
(V (hf)i V(hf) 2 ...V(hf) d ), i.e. D(hf ) := (df hf)f J=1 . Finally, we define if : (0,oo) 
(0,oo), if(r) := <p(r)/(r<p'(r)). Note that if(r) = e^C 

The following lemma provides some basic properties of the functions introduced above. 


Lemma 3.1. Let p = 1, $1 > 1, $2 > 0 and .4,(1, $ 1 , $ 2 ) be satisfied. Then, almost 
everywhere on B\, we have for the matrix-valued functions Mfff and M£ the relations 


M J 

- L VJ - a 1 


if (a )Mp + aif'(a) 


A 


VaVa T ' 
A Va T AVa A ’ 


M^° = 0 and Mfa = 0, 


and for the functions F^, F^ and F^° the relations 


F^f = if(a)F^ — a if'(a) and F^° = d — 2. 
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The relations of Lemma 3.1 can be found in [EV03], see also [MRV11], Note that our matrix¬ 
valued function M^ coincides with the matrix-valued function S^ of [MRV11], g € {tc,(r}. 
Indeed, the matrix S is dehned by 

2 {Sf)ij := -F^a lj 5ij + div(h£)a lj §ij - d k (h£)ja kl - d k (h£)ia kj + (hf) k d k a lj , 

which coincides with 2 (Mg)ij. However, since no details of the proof are given in [EV03, 
MRV11] we give a proof of Lemma 3.1 in Appendix A. 

In the following proposition we prove the quantitative estimates needed for our quantitative 
version of a Carleman estimate. 


Proposition 3.2. Let p = 1, > 1, $2 > 0 and A( l,i?i,i?2) be satisfied. Then, for all 

f 6 M. d and almost everywhere on B\ we have the estimates 

\Fp — Fp°\ < <jC' f , \i T M^\<aC M fAi, \F^\<C f , and \L 0 fi(a)\ < C^/a, 

where 

C' F := 3d^ /2 r? 2 , C F := e^ [V^i(C F + g) + \d- 2|) , 

C M := llddf /2 d 2 , Cp := (V<h{C' F + p) + d - l) . 

Proof. Recall that the coefficients aP are almost everywhere on B\ differentiable. We start 
by estimating | F^ — F^°\. By Lemma 3.1 and the definition of F^ we have 


We calculate 


F A _ F A o — 

CF <7 


cjL^g 

Vc j t AVcj 


1 - (d- 2 ). 


Va T AVa = a 2 x T A^ 1 AA Q l x and L^a = a 3 x T A 0 1 AA 0 l x — a 1 div(AA 0 1 x) 
and obtain 


F A _ F Ao = ^ div ( AA 0 ±X ) 


x T A 0 1 AA 0 1 ' 


-d = 


V 2 Ei,i e J A 0 1 x + u 2 Tr (AAq 3 ) 


x T A 0 1 AA 0 1 


x 


- d. (9) 


By our assumption on ellipticity and Lipschitz continuity we have \ejA Q 1 x\ < \fid\G, 
Ya i=i\di ail \ — and x t Aq 1 AAq 1 x > 'd^ 2 < 7 2 . Hence, 


| F* - F^° | < add\ , 2 d 2 + 


u 2 Tr(AA 0 - 3 ) 


x T A 0 1 AA 0 1 


-d 


x 


= addl /2 d 2 + \TT{T)\, 


where 


T := gAA 0 1 — Id and g : = 


<7 


x T A 0 1 AA 0 1 ~ 


x 


Next we estimate |Tr(T)| < d||T||. We have 

P II = II <7-^ A) 1 — qAqAq 1 + gAoA 0 1 — Id|| < \g\ || AA 0 1 — AqA 0 1 1| + \g — 1| 
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By 


By our ellipticity assumption we have g = a 2 / (.t t A 1 AA 0 A) — A an d 
x t Aq 1 {A 0 - A)Aq 1 x 


\9~ 1| = 


x T A 1 AA n 1 ' 


<P-A)|| 


x T A 1 A A 


1 o 

i-i 


x T A AA q ~x 


— < Ap-A||- 


i 0 x 

For the norm ||^4 q — ^.|| we use our assumption on Lipschitz continuity and obtain 


||A) — A < \/||A — AliiP — Alloo < AM < Av / A _cr - (io) 

By the same argument we have ||AAg — AqAq || < A Hence we have ||T|| < 2r? 1 / Act 

and we obtain the claimed bound | — F^°\ < <rC^. 

Now we turn to the estimate |£ t MA| < <tCm£ T A- From Lemma 3.1 we infer that 
M^A^A = 0. Note also that h^° = x and D(h^°) = I. Hence 

M* = M* - M^A^A = Si+ S2 + S3 + Sj 


where 

Si :=-\(F? - F*>)A, S 2 :=l{dw(h£oA)-dA) and S 3 := -\a {D(h^) - I) . 
For the second summand S 2 we have by the product rule for the divergence and Eq. (9) 


s A 


v 


a 2 <A 


x r AF 1 AA ~ 1 


i 0 0 


x 


T I 

AA^x + A{Ff - F*>) ] =: 5 S + S 6 . 

- i,j =1 


Hence, Si and A cancel out and we have = 5s + S 3 + S’J. For S 5 we calculate 

S 5 = \s 7 + ^S 8 + Sg, 

where 


A := 


er 


x T A 0 1 AA 0 
a 2 A 


31- [(Va^A^x]^ 

n •L 


A := A TrHri )2 Z x ’ Tj4 o ’(AAA A e feM) A, 
(*■ A ^A x ) k 


Sg := A — 


a 


2 a 


(x^A^AA^xF TA ° 1AA °' AA ° 1x ' 


Here dkA denotes the matrix (AcA)A = 1 . For the Frobenius norm of S 7 we have using 
g < ft 2 , Y.j\dka lj \ 2 < 'd ‘ 2 and \AA^ l x\ < i9^ 2 a 


IIAIIf = g y 


Y I (VaA T M 

\ i,j =1 


0 A| 2 < 5a 


W 


Z |AA| 2 |(HA A)fc | 2 < Vdd^ 2 d 2 a, 


k =1 
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and hence we obtain 


|£ t S 7 £I < \\S 7 U t Z < ||S 7 || F £ T £ < ||5 7 ||F^ T Ae < Vd-&l /2 -&2<rfM- 

For Sg we use Holder’s inequality and |HHq 1 x| < y/Wix' 1 Aq 1 AAq 1 x and obtain 

\ 1/2 


ie T 5 8 ei < , (E^^A)^ 


-^l 2 


(,t t h 0 1 aa 0 1 x ) 3 / 2 

Since ||<9fcH|| < V / d||9fc^4||oo < Vd '&2 and x^ A^ 1 AA^ l x > ■d]” 1 x T Hg 1 Aq l x we obtain 

|£ T S 8 £| <- f dd 2 x T A- ] A- 1 x < dJ 2 < di?i /2 7? 2 ^ T ^- 

1 ^ - (x^A^AA^xf/ 2 0 0 ~ (x t Ao 1 Aq 1 x ) 1 / 2 1 

For Sg we calculate 

licTx)T{v - x) - {s ~ xfCxl K - 

where we used the notations C = Aq 1 AAq 1 xx t Hq 1 and y = AAq 1 x. The inequali¬ 
ties \Cx\ < cr 2 'd^ //2 (x T HQ 1 HHg 1 x) 1 / 2 , |C' T x| < cr' 2, d^ 2 {x T A^ 1 AA^x) 1 / 2 , \x — y\ < ||H — 
Ho||(x T Hg 1 Hq 1 x) 1 / 2 , and x^A^AA^x > i?^ 1 x T Hg 1 Aq l x imply 


\fS 9 £\ < 


2e T ^u 2 ^ /2 p-H 0 ||(x T H 0 - 1 H 0 - 1 x) 1 /2 < 2 eA^Hl\\A-A 0 


(x T H 0 1 AA 0 1 x ) 3 / 2 
Using x t Hq 1 Hq 1 x > '( 9 ^ 1 ( 7 2 and Ineq. (10) we obtain 

\H T Sgti\ < 2fM^\\A - Ao\\ < Z^A^ 2 ^* 




For £5 we finally obtain the estimate 

|£ T S 5 £| < 3 d^ /2 ^ T AC. (11) 

Now we start estimating 53 . We have D(h£) = gD(AA^ 1 x) + diag(V<?)XHg 1 A, where 
X := (x ■ ■ ■ x) T and diag(v) := diag(r»i,..., va) for v £ M rf . For dkg we calculate 


dk9 


^ 0 lx 

,,T 4-1 4 4 -I™ 

<X/ VIQ VT.VJ.Q Jb 


A ^ T 

Jb VT.Q Jb 

(x T Aq 1 AAq 1 x) 2 


(2ejA 0 1 AA 0 1 x + x T H 0 1 (d k A)A 0 1 x) . 


Hence, 

S 3 = ~A(S W + S ll + S 12 ) 

where 


5 io := 9 D(AA^x) - I, 

s n ■= T ._ 1 , ._i [A 1 ® - gA- ] AA- l x] x T A^A, 

X vj.q vAvAq 
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*-’12 :=- T A-lAA-l (^A) 1 ( d i A ) A 0 1 *)?= 1* T A) lA 

aHA, HA x 


i 0 

For the norm of ,S'io we calculate using D(x) = / 

i-i 


||S 10 || < g\\D(AA^x - ®)|| + \g - 1| = g ||A„ :l A + C-l\\ + \g-l\ 

< sll A 1 ll ■ 11^- — All + g lie’ll + \g — 1|) 


where 


c = {±W*)eW*) ij=i . 


k =1 


The existing bounds g < i9 2 , \g — 1| < ||H — Ho||, ||H — Ho|| 5 AvA 0- and 


lie’ll < IICIIf < 


N 


d 

y: |<9jaJ fe | 2 |e^o x x| 2 < 




1 


E^ 1 *! 2 


< y/dzd\d 2<T 


give us ||S'io|| < 3v / di?y 2 'd2 cr - For 5n we calculate using Cauchy-Schwarz’s inequality, 
\AAq l x\ < \/?9i(x T j4q 1 AAq 1 x) 1/ ’ 2 , x t Hq AHq A > Aq 1 x, T = gAA^ 1 — I, and 

||T|| < 2 '& 7 1 ^ 2 '&2^ as before 


ll-S'iiH < 


2 VW^ 

(x T Ag 1 AAq 1 x ) 1 / 2 


\{I-gA^A)A^x\ < 20 1 ||T|| < 4^ /2 i? 2 a. 


For S \2 we use Cauchy Schwarz’s inequality and |^4A 0 1 x| < \Adi(x T A 0 1 AA 0 ' 1 i) 1 / 2 , and 
obtain 


IIS12II < 


y/tijg 

(x T Ag 1 AAq 1 x ) 1 / 2 


| (x T ^ 0 1 (a^)^ 0 1 x)f =1 


Since |x t Hq A A lx l < -v/dAa^A lx we Fave |(x t HA (<9iA A 1;c )f=il < cM 2 % T 
A^A^x. Hence, using x T A^AA^x > d^x 1 - Aq 1 A^ l x, g < d 2 and x t Hq 1 Hq A < i?icj 2 
we obtain 

||5i 2 || < digdd 2 {x T A^ x Aq 1 x) 1 / 2 < dd 7 / 2 d 2 cr. 

Note that for any matrix C we have by ellipticity of A the estimate S^ ACS < i?i||C'||£ t .A£- 
Hence, we obtain 


lA^I < ^ T AA||Ao + 5 n + 5*1211 < 4dd| 1/2 i?2 afA£. (12) 

From Ineq. (11) and (12) we conclude that \S T M^S\ = |C T (5s + S 3 + Sj)£| < < 7 CmS T AS- 
In order to prove the bound \F^\ < Cf we use the formulas = ip[a)F^ — aip'(a) and 
F^° = d — 2 from Lemma 3.1, our previous bound | F^ — F^°\ < aC' F , ip'{ a ) = gtp(a), 
ip (a) < e^^F and a < y/d 1 on Hi to obtain 

\F*\ < 4(a) (| F? - | + | F*> - ag\) < C F . 
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For our last bound |LoV’(°')l < Cy, jo we note that V^cr) = pe^a 1 A 0 1 x, use the product 
rule for the divergence and Eq. (9) to obtain 

, , , ( (\ T , ,, . A 


4. Proof of Theorem 2.1 

As mentioned in the introduction, the proof of Theorem 2.1 is based on ideas from [EV03, 
BK05]. One might think that once the quantitative bounds from Section 3 are at hand, the 
result is obtained simply by following the existing proofs [EV03, BK05]. However, due to the 
fact that we want to have an explicit dependence on all the parameters and also treat the 
case k = 1 with additional gradient term in the lower bound at the same time (in the sense of 
Ineq. (1)), we need to employ non-trivial extensions of the techniques used in the mentioned 
papers. More precisely, the proof of Carleman estimate in [EV03] has the gradient term in 
the lower bound which is sometimes needed in the applications, but the non-scaled version is 
valid on B K \ {0} with undetermined k < 1, and the proof does not allow for a scaling which 
keeps the constants explicit. The Carleman estimate in [BK05], on the other hand, is valid 
on the whole punctured unit ball in the non-scaled version, but the gradient cannot be easily 
preserved without breaking the proof. Our extension of their techniques addresses this issue. 

First we prove a special case of the theorem and assume that 


u G C^°(B i \ {0}), p = 1, b, c = 0, and u is real-valued. 


(13) 


Note that in this case we have L = Lq. The general case then follows by regularization, 
scaling, the fact that Carleman estimates are stable under first order perturbations, and by 
adding the obtained Carleman estimate for Re u and Im u. These steps are carried out in 
detail at the end of this section. 

We set / = w~ a u. The following lemma can be found in [EV03, MRV11], For completeness 
we give a proof in Appendix B. 

Lemma 4.1. For all a > 0 we have 



(14) 


We start the proof by providing a lower bound on the first term of the right hand side of 
Ineq. (14). With the notation 



AVf = V/ 


Vw T AVw AVf 
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we have 


/ VcrVcr T 
\ Vcr T AVcr 



T 

AVf = 0, 


and using the first identity of Lemma 3.1 we obtain 


V/ t M^V/ = iP(a)Vf T M?Vf + a^'(a)Vf T AVf 

= ip(a)V f T f + a^(a) (vf + ^ ^ AV/) AWf 

= iP(a)Vf T M?Vf + atp' (cr)'V f T AV f . 

Now using M^'S/cr = 0 established in Lemma 3.1 and the fact that is symmetric, we 
obtain 

V/ t M^V/ = 'ip(a)V f T f + a^(a)Vf T AVf. 

From Proposition 3.2 we obtain for almost all x G B\ 

V/ t M^V/ > (a^(a) - V>( o)oCm ) V/ T 7lV/. (15) 


Using Ineq. (15) we provide a lower bound on the first two terms of the right hand side of 
Ineq. (14). Lemma 3.1 implies 

F^ = i,(a)(d-2) + i;(a)B^-^(a)a, where B?:=F?-F*>. (16) 

We use Ineq. (15), Eq. (16), L 0 (/ 2 ) = 2fL 0 f - 2V/ t 7LV/, 


v/ t av/ 


Vf T AVf + 


(\7w T A\7f) 2 

Vw T A\7w 


and Green’s theorem, i.e. f uLqv = f Vu T AVv for u,v G C 2 (B\), to obtain 


(17) 


h ■= 4 a J (V/) t M^V/ -a I F^L 0 (f) 

>4 aj (aip'(a) “ i>{(y)aC M ) V/ T AV/ - (d - 2) a j f 2 L 0 ip(cr) 

-2a j (, mB. J (/W - V/ T W/ - 

We apply Eq. (34) to obtain 


h > 2a J (crip'(a) - 2i’(a)aC M + ip(a)B^) \7f T AVf 

- 2a 3 J (iI>(<t)b£ - ip'{a)a) Vw f 2 - R u (18) 

where 

Ri := (d-2)a J f 2 L 0 ip(a) + 4a 2 J (ip(a)B^ - ip'(a)a) fF> f VW J^ VW 
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-2 a J (ip(a)Bff - 'ip'(a)a) 

From Lemma 4.1 and Eq. (18) we infer that 


( Vw T A\7ff 


Vw T A\/w 


- w~ Q fL 0 u 


Ri + J y u {r A Tj w ( w a L 0 u) 2 >2a J (aip'(a) - 2ip(a)aC M + ip(a)Bpf) Vf T AVf 

\ / w T A\ / w 


- 2a 3 [ (ip(a)B^ - ip'(cr)cr) Vw + 4q ,2 




te- 


xxr 


We use the identity ip’ = /j,ip and the estimate \B^ \ < Cm& from Proposition 3.2 and obtain 
with the notation C^ := /x — 3 Cm 


R\ 


+ J 2 2 “ C * J <^MV/ T W/ 

+ 2^cJ^(a)^A^f + ^ 


X7w t AX7w 2 


x/r 


£>/■ (19) 


Note that /x — B^/a > /j, — 3 Cm = CL > 0 by our assumption on /x and Proposition 3.2. The 
identities (17), 

Vf T A\7f = w~ 2a 'Vu T AVu + a 2 w~ 2 f 2 Vw T AVw - 2aw~ a ~ 1 fVw T AVu 
= xW 2q Vxx T j4Vxx — a 2 w~ 2 / 2 Vxx; T AVxx; — 2aw~ 1 fV w T A\7 f 


and 


imply 


fVw T AVf = fw~ 1 D f Vw T AVw - ^f 2 w~ l F^Vw T AVw 


Vf T AVf = ur 2a Vu T AVu - a 2 w~ 2 f 2 Vw T AVw - 2afw- 2 D f Vw T AVw 


+ axe 2 f 2 F^\7w T AVw — 


(\7w T A\7 f) 2 
\7w T AVw 


Combining this with Ineq. (19), and using a > w, ip > 1 and the bound Ffj > —Cf from 
Proposition 3.2, we derive at 


w 


R2 + J - 2aC » J w ~ Za+ ^AVu 

'Vw T A'Vw 


— 2C^CfoT / aip(a)- 


w“ 


f2 + ia2 f y»^yu, D?i (20) 


xxr 


where 


R 2 ■■= (d-2)a J f 2 L 0 ip(a) + 4a 2 J (ip(a)B^ - ip'(a)a + C /J ,aip(a)) fD f 


Vw t AVw 


w~ 


-2a (ip(a)B^-ip'(a)a) 


(Vw T AVf) 2 

X7w T AVw 


- w a fL 0 u 


fa ,(Vxe T AV/) 2 
+ 2 C>/w>M Vw r AVw ■ 
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We drop the (positive) term 2ctC^J cn/?((T)V/ T HV/ in Ineq. (19) and add the inequality 
obtained in this way with Ineq. (20). This gives us 


R 


Vw r A\/w 2 


Dj, (21) 


+ / - aC » / “ _2 ° +lv >‘ T - 4v “ 

+ a 2 (a - + 4a 2 

where R := (R\ + R 2 )/2 is given by 

R = (d-2)a j f 2 Loi/j(cr) + 4a 2 J (ip{a)B^ - ip\a)a + ^aip(a)^j fD f ^ W 

-2a J (i,(a)BS - $(„)a - ^IP/w 

+ 2 a J (ip(a)B^ - ip'(a)a) w~ a fL 0 u. 

Next we provide an upper bound on \R\. Inequality \B^ \ < Cm& from Proposition 3.2 and 
our assumption /j, > 3 Cm imply the estimates 



— 

(7 2 



— 




B£ _ Ip _ 3 Cm 
a 2 2 

~~ — ^A 4 + \c M 


4 

- 3 ^’ 

. 3 l n . 3 


and | B^/a — y\ < 4/i/3. Hence, using \Loip(a)\ < C^/a from Proposition 3.2, w < a, 
ip' = nip, a < ni w i V’( (J ) < and Vrc T ylVtt! < (ip'(a)) 2, d\ < d 2 we obtain 


|i?| < Cibda /-h 

./ w 


< I\ 


P 

w 

f 2 


J + + » j ' “l/ L o“l 


« / - +a » /IhH + a +a 


u; 


u> 


Vw T j4Vw 


where 

K := max jdC^, 6^ie M ^U? 2 j < Qd/ini^'^d\ (\f®i(C' F + n) + d^j . 

Using (a — 6) 2 < 2a 2 + 26 2 , S/w^AS/w < df and w < a < y/d\ on B\, we obtain 

w (Vv/^AV/y (Vw T W/) 2 p} , eg./ 2 / 5 -U+ 

7 2 Vu> T HVu; (VuTHVu;) 2 re 2 to 1 w 2 2 w ’ 

Furthermore, for all t > 0 we have 

a J w l ~ a \f Lqu\ < V ^" a J PP + \ J w 2 ~ 2 a \ L o u \ 2 > and 
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W 1 \fDj\ < t J w V 2 + J W 2 D 


- 2 n 2 

/• 


Hence, 


\R\ 


ot ( 1 + cat + 


Cp'&i aty/&i\ 


w - 1 / 2 


F)/ 

+ a ( a ^l+2^ 2 ) j w^D} + \fw^\L„u\ 2 - (22) 

From Ineq. (21), cn/’( fJ ) > w, the bound Vw T AVw > $]“ 2 /hj” 2 e _2/iV ^i and Ineq. (22), we 
finally obtain 

K\ J w 2 ~ 2 a \L 0 u \ 2 > K -2 J vj~ x f 2 + K 3 J w~ 2 D) + IU J w~ 2 a+ 1 Vu T AVu, (23) 


where 


I<i := K/2 + (i?i^i) 2 e 2 ^, 


k 2 := a 


p-2/iv^r 

e o 

At ^ a - 


(tfl/il) 2 


,-2/VA 


/Vi : = 


4--Ji 

(tVi) 2 


o ; 2 — 2 K'&'y ' 2 a and 


cr+- S 
+ 

' i^l 

o 2 - K 

\ c F$i 1 
1 + ^ 1 

tc 

1_ 


2 


a, 


At 


K4, := Cfj,a. 

Now we choose t large enough, such that the coefficient of a 2 in A 3 is positive, and thereafter 
we choose a sufficiently large, such that K 3 is non-negative. Our particular choice is 

t = h:= Ix/f 2 and a > K'd’f 2 =: ai . 

8 


Then 


K 3 = 2 


s -2 


r a' - 2Kiff 2 a > 0 


and Ineq. (23) implies 

K\ j w 2 ~ 2 a \L 0 u \ 2 > K 2 J w~ l f 2 + A 4 J w~ 2 a+ 1 Vu T AVu. 

If additionally 


(24) 


f/ , / q 2 2 r _ 

o > — V \ 3 4-- —: 021 

p U jr p 


where 


g — 2fJ.y/Wl 

p := C “ (W 


j Q ■= CuCf 


^-2iisf&x 

('l?l^l) 2 


+ A ( t\ + 


\f &1 


r := A 
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then 


1 -2/V$7 

A ' 2 * 2 C " (W ^ =: A ' S “ 3 ' 

Moreover, it can be shown that aq < « 2 - By using min {K^,Cn} = we obtain from 
Ineq. (24) for all a > bo 

ur 2 a+ 1 Vu T AVu < C j w 2 ~ 2 a (L 0 u ) 2 , (25) 

where 

do = &o(d, , $ 2 , /i) := maxjai, 02 } = a 2 and C = C{d,d 1 ,^ 2 , p) ■= ~r2- 

B- 5 

This proves Theorem 2.1 in case (13). 

Now we lift the restriction that u G C“(5i \ {0}) and assume that u G bF 2,2 (M d ) real¬ 
valued with support in B\ \ {0}. Let (f> be a non-negative and real-valued function in C£°(M. d ) 
with the properties that ||</>||i = 1 and supp^> C B\. For e > 0 we define (f) £ : M rf —> Mg" by 
4> e (x) = £~ d <f>(x/e). The function cj) e belongs to and satisfies supp^ e C B e . Now 

define 

u £ = (j) £ *u, u e (x)= <j) £ {x - y)u(y)dy. 

J R d 

If e is small enough, say £ G (0,£o), then u £ G C%°(Bi \ {0}), see [Zie89, Theorem 1.6.1]. 
Hence we can apply Ineq. (25) to the function u £ . By definition we have 

d d d 

L 0 u e = -^2 di{a lJ dj((l) £ * u )) = -^2 ( d i a%3 ) (</>e * djU ) - ^ a 13 ( (j) £ * didjU ). 

*,l=i *,l=i *,l=i 



Since djU and didjU are elements of L 2 (M d ), [Zie89, Theorem 1.6.1 (iii)] tells us that (j) £ *djU — > 
djU and cj) £ * didjU —> didju in L 2 (M rf ) as e tends to zero. Hence Lqu £ —> Lqu in L 2 as £ —> 0. 
Since obviously Vu e —> Vu and u £ —> u in L 2 as £ —> 0, and ie _1 is bounded both from below 
and above on suppu £ uniformly for all £ G (0,£ 0 ), we obtain Ineq. (25) for u G IT 2 ’ 2 (M rf ) 
real-valued with support in B\ \ {0} in the case b, c = 0. 

Now we lift the restriction on p and assume that p > 0 arbitrary and u G VF 2 , 2 (M <i ) real¬ 
valued with support in B p \{ 0}. We introduce the scaled coefficient functions a 13 (x) := a lJ (px) 
on B\, and define the scaled elliptic operator 


Lq : — 


d 


*,l=i 


Obviously, the Lipschitz constant of Lq is p $2 and the ellipticity constant is i?i. By our 
assumption on p, we can apply Ineq. (25) to the function u : B\ — > M, u{x) := u(px), and 
obtain for all a > do 

a J w 1 ~ 2 a Vu r A'Vu +a 3 J w~ l ^ 2 a u 2 <C J w 2 ~ 2 a (Lou) 2 , 
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where 


a 0 := a 0 (d,i}i,p'&2, p), C := C'(d,i?i, p&i, p) (26) 

and w(x) := w(px ) for x E Ri. Note w = w P)fl = ip(a(x/p)) and hence w(x) = wi^(x). By 
the change of variables y = xp, we obtain for all a > do 

ap 2 j u; 1 _ 2 a Vii T AVu + a 3 j w~ l ~ 2a u 2 < Cp 4 J w 2 ~ 2a (Lqu) 2 , (27) 

which proves the theorem in the case u E B / 2 , 2 (M <i ) real-valued with support in B p \ {0} and 
b, c = 0 . 

Now we lift the restriction that u is real-valued and, i.e. we allow complex-valued functions. 
We apply Ineq. (27) to the real part Reu and imaginary part Imtt and add these two 
inequalities. We obtain Ineq. (27) for u E IR 2 , 2 (M rf ) complex-valued, since V and Lq commute 
with Re and Im, A is positive, and (VRe«) T d(VReM) + (V Im u) T cl(V Im u) = Vu T 7lVu. 

Finally we lift the restriction on b and c. Let u E IF 2 , 2 (M rf ) be complex-valued with 
support in B p \ {0} and let 6 , c E L°°(B p ) be arbitrary complex valued functions. From 
|ai + ci 2 + as\ 2 < 3(|ai | 2 + |a. 2 1 2 + 1 0 , 3 1 2 ) for complex numbers a* we obtain 

\Lqu \ 2 < 3 (\Lu\ 2 + | 6 T Vtt | 2 + |eu| 2 ) < 3 (\Lu\ 2 + WbW^PiVu 1 AVu + HcH^Jrtl 2 ) . 

By Ineq. (27) and w{x) < 1 on B p , it follows that for all a > do we have 

l ~ 2a \u\ 2 < 3Cp 4 J w 2 ~ 2a \Lu\ 2 

Now we can subsume the two last terms on the right hand side into the left hand side by 
choosing a sufficiently large. In particular, setting 

C:= 6C and a 0 := max |a 0 , Cp 2 \\b\\ 2 00 ‘&\^, C 1 ' / 3 / o 4 // 3 ||c||^ 3 v / ^i| , (28) 

we obtain the statement of the theorem for all a > ao and all u E IR 2 , 2 (M rf ) with support in 

B P \{ 0}. □ 



A. Proof of Lemma 3.1 

First we prove the fourth relation. We use Vw T AVw = cr~ 2 p'(a) 2 x T Aq 1 AAq l x. 


- Lqw = p"(a)-^x T A 0 1 AA q l x - p'{a) 


■■ t A 0 l AA 0 l x div (AAq l x ) 


( 7 ° 


o 


and since p = 1 


and obtain 


, w wtp"(a) 

oi\) (cr) = 1-——-. 

ap'(a) p (cr) 2 


F A = -1 w <P 

x in x 1 


'(cr) w wa div [AA 0 1 x) 


ip'(a) 2 ap'(a) ' p'(a)x T A 0 X AA, 1 ™’ 


(29) 


(30) 

(31) 
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Note that w = <po a and Eq. (29) with <p replaced by the identity reads 


— Lo<r = —a 3 x T A 0 1 AA 0 1 x + a 1 div ( AA 0 1 x) . (32) 

Using Va T AVa = (J ~ 2 x T A ^ 1 AA ^ 1 x and Eq. (32) we have by the definition of 

dW(AA^ l x) = xT t 4 o^V^ (i? A + 2 ). (33) 

From Eq. (30), Eq. (31), Eq. (33) and the definition of iji we obtain F = — ai/i'(a) + i/i(a)Fi 1 
i.e. the fourth estimate of Lemma 3.1. 

The fifth identity of Lemma 3.1 F^° = d — 2 is a direct consequence of Vcr T AVcr = 
(t~ 2 x t Aq 1 AAq 1 x, cr 2 = x t Aq l x and Eq. (32). 

We now turn to the proof of the first relation. By definition we have Mi = S 1 + S 2 + S 3 + &J 
where 

Si ■= ~f£A, S 2 ■■= ^ div (h£ o A) , and S 3 := ~AD(h£). 

For the first summand we have, using the fourth identity of Lemma 3.1, 

Si = -1 {iI>(<t)F* - ai/i'(a)) A. 

For the second term we calculate using the chain rule and the product rule for the divergence 



/ p(cr) aAVa 
\ crip'(a) Vcx T AVer 


-o"il>'(a)A + -ip(cr) div 


aAS/a 

Vcr T AVcr 



For the third summand we use the chain rule to see that D(h (j) = D(i/i(a)hi). Hence, by 
the product rule 

S 3 = ~il>(<r)AD(h£) - ^Ai/i'(a)(Va)(hi) T . 

Putting everything together we obtain that 


Mi 


+ cnl)'(cr)A 


+ cnli'{a) 


A 


fp-(A(Va)(h^) T + h^(Va) T A) 

1 ( „ , VctVcj- 1 - A 

- - A————A + A — . — A 

2 \ VcAAVct Vcr-'-AVc ) 


To prove the second relation M^° = 0 we have by h£° = x 

div (x o Ao) = dAo and D(hi) = I. 

Hence, using F^° = d — 2 we obtain 

= -\f^A 0 + 1 div o A 0 ) -\A 0 D(hi) - lD(hi) T A 0 = 0. 

For the proof of the third relation we infer from the proof of Proposition 3.2 that = 
(M^- — Mi 0 Aq 1 A) = (£5 + £3 + Sj). Moreover, we have that 2 £5 = £7 + Sg + 2£g and 
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2S% = —^4(5io + <511 + S' 12 ). We rearrange the terms, recall that g = cr 2 /x T A 0 1 AA 0 l x and 
dkA = (<9fc(W)f J=1 , set v = (x T Aq l (dk.A)AQ 1 x)^ =1 , use D(AAq 1 x) = A^ X A + C, where 


C 


d 

(j^^a^kA^x) 


d 


1 

i ,]=1 


and obtain that M^ = T % with 


Ti ■= § [(VaV) T AA^x]l j:ml - |C t A 


9 n T 


To := ^ 9 f - ^vx t AT 1 A - lAC, 

2 x t Aq 1 AAu 1 x 0 2 

j 1 .. —_^_ t 1 A~ l A A~ 1 A A~ 1 r A _ - _4 4 _1 tt T /4” 1 44 _1 /1 

3 ■“ rTA-iAA-K 0 0 0 + t t 0 0 0 ’ 

ily -TIq 7T.X1lq i/y <ly YiQ JlJl-Q Ju 


9 


Ta T A-i a a- 1 ^^0 1 AA 0 1 xx T A 0 1 A gAA 0 L A, 

X 1 Aq ^4y4g x 


Ta := 2 A - 


n ■= 


—r- 1 — AAr. 1 xx T A n 1 A, 

x^A^AA^x 0 0 


9 


2x T A 0 1 AA 0 1 ; 


-y4A n 1 xf T y4 — 


5-4 


2x T y4g MAq "X 




, v jaja n x. 

-1 A A-l„ U 


If we multiply from the right with Va = a 1 Aq 1 x , then one easily sees that Tj = 0, 
i 6 {3,4, 5, 6}. It remains to show that T{S7a = ToSo = 0. For this purpose, we calculate 


{Vtf^AA^x 


- i,k= 1 


An 1 ' 


X = 


dia jk eJAA 0 1 xe%A 0 h 


A k i 

= AA^xCij = (C t AAq 1 x)j, 


hence [(Va^^A^g 1 x] 1 * A 0 l x = C 1 -AA Q l x and T\\7 a = 0. 

For T 2 V a we calculate 

d 

(CAq 1 x)j = ^2 ejA q 1 xdiai k ej: A^ 1 x = x T y4g 1 (diA)AQ l x = 

j,k =1 

and obtain 

r 2 Vu = - ^-ACAn L x = J!~A(v - CAn l x) = 0. □ 

2(7 2(7 2(7 

B. Proof of Lemma 4.1 

We recall that / = w~ a u and calculate 

—w~ a L()U = a(a — l)w~ 2 fV w T AVw + 2cra; -1 V/ T y4Vu; — afw -1 Lqw — LqJ 
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(34) 


Hence, 


where 


h = 


= —Lof + a 2 fw 2 Vu> T AVw + 2aw 2 (Vu> T A\7w)Df 
n Vw T AVw _ 9 


4a 


P : = 


-D f — AaDfLof + 4a 
w z J w z 


3 VwUlV» +f . 


(35) 


2 Ww t AVw r 

a 2 -/- 

w 


w 


Ww T AX7w 


Lof) > 0 . 


We use Green’s formula, i.e. f uL$v = f Vu T AVv for functions u,v £ Cq(B\), and obtain 
for the third term in Eq. (35) 


\/w t A\/w 


/ V7,.. T/tvyr f f 2 r „„ f f Y 7 „„T 




D,f- I I 


w 


f W,„T 


Vw l AVw 

2 \T 


w 


VroMV/- 


/ 2 u )2 VtyT ^ Vu; " / V (£) ^ Vw - 


By the quotient rule, the last term equals to the sum of the two first terms. Hence, 

f X 7 w t AX 7 w 


ur 


Dff = 0. 


For the second term in Eq. (35) we have 


4a j D f Lof = Aa j + 2a j F wf L of- 

By Green’s formula, we have 

2a J F^fL 0 f = 2 aj V(i^/) T HV/ = 2 aj /V/ T HVi^ + 2 aj i^V/ T HV/ 
= aj (Vi^) T HV(/ 2 ) + 2 aj i^V/ T HV/ 

= aj i^L 0 (/ 2 ) + 2a J i^V/ T HV/. 

Hence, the second term in Eq. (35) reads 

Vw T AVf 


4aJ D f Lof = 4c,Jw VwTAVw 
From Rellich’s identity [Necl2, Eq. (5.2)] we obtain 


L 0 / + 2a / i^V/ T HV/ + a / F^L 0 (f z ). (36) 


J (/4) t V/L 0 / = ~J V/ t PV/, 


(37) 


where 


h* = 


wAVw 

Vw T AVw 


and B := di v(h^A) - AD(h£) — D(h^,) T A. 
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From Eq. (36) and (37) we obtain 


4a J D f L 0 f = -2a J Vf T BVf + 2a j F^Vf T AVf + a j F*L 0 (f 2 ). 
By definition of this gives 


-4a J D f Lof = 4 af (V/) T M^V/ -aj F*L 0 (f 2 ). 


□ 
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